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Foreword 


This  paper  was  prepared  for  presentation  at  the  SNOW  BALL  Synrposium 
held  at  Sandia  Base,  New  Mexico,  on  19-21  JanuEury  19^5 >  under  sponsorship 
of  the  Defense  Atoicic  Support  Agency  (DASA)c 

The  work  discussed  was  sponsored  hy  DASA  as  part  of  the  Nuclear 
Weapons  Effects  Research  (NWER)  program  Subtask  I3.OO9,  Dynamic  Res-ronse 
of  Soils.  Additional  funding  support  was  provided  by  the  Office,  Chief 
of  Engineers,  U.  S.  Amy,  as  part  of  the  Military  Engineering  and  NWER 
program. 

The  studies  described  were  conducted  under  the  general  direction  of 
Mr.  F.  R.  Brown,  former  Chief  of  the  Nuclear  Weapons  Effects  Division; 

Mr.  G.  L.  Arbuthnot,  Jr.,  Acting  Ctiief  cf  the  Nuclear  Weapons  Effects 
Division;  and  Mr.  L.  F.  Ingram,  Chief  of  the  Physical  Sciences  Branch. 
During  the  preparation  of  this  report  Col.  John  R.  Oswalt,  Jr.,  CE,  was 
the  Director,  and  Mr.  J.  B.  Tiffany  the  Technical  Director  of  the  Waterways 
Experiment  Station. 

The  author  wishes  to  express  his  appref'iation  to  Mr.  Paul  Hadala 
for  his  Invaluable  contributions  toward  the  successful  conpletlon  of  the 
wave  propagation  test  conducted  on  an  undisturbed  sanple  of  soil  from  the 
Suffield  Experimental  Station  in  the  Waterways  Experiment  Station  Sknall 
Blast  Load  Generator.  The  author  ia  also  indebted  to  Mr.  D.  W.  Murrell 
and  Mr.  J.  K.  Ingram  for  their  cooperation  in  making  a^i-ailable  ground- 
motion  measurements  without  which  this  paper  could  not  have  been  written. 
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Summary 


In  this  paper  the  dynamic  properties  of  the  soil  profile  near  the 
location  of  the  Snowball  detonation  (Suffield  Experimental  Station, 
Watching  Hill  Site,  Alberta,  Canada)  are  given  which  pertain  to  the 
prediction  of  the  airblast-induced  ground  motions  in  the  superseismic 
region.  Predictions  based  upon  the  available  dynamic  properties  are 
given  and  compared  with  the  ground  motions  measitred  by  the  Waterways 
Experiment  Station. 


CORRELATION  OF  OPERATION  SNOWBALL  GROUND  MOTIONS 


WITH  DYNAMIC  PROPERTIES  OF  TEST  SITE  SOILS 


Introduction 


1.  The  purposes  of  this  paper  are  to  present  the  results  of  dynamic 
tests  which  were  conducted  on  undisturbed  samples  of  soil  obtained  near 
the  location  of  the  Snowball  shot,  and  to  show  how  the  soil  properties 
were  used  to  predict  peak  transient  displacements  In  the  superselsmlc 
region.  A  comparison  of  the  coinputed  and  observed  ground  displacements 
at  stations  200  and  250  ft  from  ground  zero  (GZ)  is  given. 

Site  Conditions 


2.  The  Snowball  shot  was  detonated  at  a  location  on  the  Suffleld 
Experimental  Station,  Alberta,  Canada,  known  as  Watching  Hill.  Within  the 
Immediate  area  of  the  explosion  the  ground  surface  was  essentially  level 
at  an  elevation  of  approximately  216?. 0  ft  (mean  sea  level).  A  general 
description  of  the  geology  Is  given  In  reference  1.  The  location  of  GZ 
was  In  the  southern  end  of  an  area  known  as  the  Ross  Depression.  This 
depression  has  apparently  once  been  covered  by  a  large  lake,  and  the  soils 
to  a  depth  of  200  ft  are  lacustrine  deposits  consisting  of  fairly  uniform 
beds  of  clays  and  silts  with  occasionail  sand  lenses.  As  is  usually  the 
case,  desiccation  and  weathering  have  altered  the  upper  30  ft.  A  recent 
seismic  survey  indicates  that  the  upper  23  ft  of  this  deposit  has  a 
seismic  velocity  of  pproxlmately  1050  fps,  whereas  from  a  depth  of  23  ft 
to  200  ft,  references  1  and  2  indicate  a  seismic  velocity  of  approximately 
5100  to  5^00  fps.  The  water  table  was  located  at  a  depth  of  23  ft  at  the 
time  of  the  test. 

3.  The  bedrock  at  the  site  is  also  described  by  C.  H.  S.  Joms,^  but 
in  general  the  surface  of  the  bedrock  at  a  depth  of  200  ft  is  too  deep  to 
Influence  the  ground  motion  in  the  superselsmlc  region  which,  of  course, 

is  the  subject  of  this  study. 


*  Raised  muaerals  refer  to  items  in  the  Literature  Cited. 


Soil  Profile 


Exploratory  borings 

k.  The  subsurface  conditions  of  the  site  were  delineated  by  borings. 
The  majority  of  these  borings  were  used  for  construction  of  sand  columns 
and  installation  of  free-field  transducers,  llie  locations  of  borings  in 
which  undisturb  3d  samples  were  taken  for  i  sting  by  the  Waterways  Experi¬ 
ment  Station  (WES)  are  shown  in  f^g.  1  and  are  designated  borings  A,  B,  C, 
D,  E,  and  F.  WES  also  obtained  undisturbed  samples  from  the  borings 
designated  R  and  N17  for  the  Air  Force  Weapons  Laboratory;  the  results  of 
dynamic  tests  and  classification  tests  conducted  on  these  sasqples  have  been 
reported  by  Davisson  and  Maynard."^  Only  one  continuous  undisturbed  boring 
(boring  A)  was  drilled;  a  relatively  few  undisturbed  samples  were  obtained 
at  various  d^ths  from  the  other  borings. 

Laboratory  test  results 

5.  The  samples  sent  to  the  laboratory  were  subjected  to  routine 
identification  and  classification.  Unconfined  conipression  strengths  were 
determined  and  several  consolidation  tests  were  conducted.  Vibration  tests 
and  unconsolidated,  undrained  triaxial  tests  were  conducted  and  reported 
by  the  firm  of  Shannon  and  Wilson.  In  addition  to  the  tests  described 
above,  a  new  type  of  dynamic  test  was  conducted  on  one  of  the  Suffleld 
samples  in  the  Snail  Blast  Load  Generator  (SBLG)  at  WES;  this  test  will 

be  described  in  detail  later. 

6.  Fig.  2  shows  typical  grain-size  distribution  curves  obtained 
frcan  representative  samples  from  various  depths.  As  a  result  of  the  vari¬ 
ous  laboratory  tests  described  above,  the  soil  profile  has  been  divided 
into  representative  layers  as  shown  in  fig.  3-  From  the  surface  to  a  depth 
of  approximately  12  ft,  the  soil  is  a  friable  silty  clay  with  a  water 
content  of  approximately  9  percent  and  an  unconfined  cocpresslve  strength, 

,  ranging  from  3*75  to  4.5  tons/ft  .  Atterberg  limits  determined  for 
saaples  from  this  stratum  indicate  that  the  liquid  limit,  LL  ,  ranges  from 
38  to  43  and  the  plastic  limit,  PL  ,  from  21  to  23;  therefore,  this  soil 
would  be  classified  as  CL  in  the  Unified  Soil  Classification  System. 

7.  From  a  depth  of  12  to  27  ft  the  soil  is  a  brown  silty  clay  with 
a  water  content  of  35  percent  and  an  unconfined  coaga*essive  strength 
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2  2 
varying  from  1.8  tons/ft  near  the  top  of  the  layer  down  to  1  ton/ft  at 

the  bottom  of  the  layer.  Ihe  liquid  limit  of  this  layer  ranges  ftrom 

48  to  52  percent  and  the  plastic  limit  is  about  21  percent.  Therefore, 

this  soil  is  somewhat  more  plastic  than  the  layer  immediately  above,  and 

it  is  borderline  between  a  CL  eind  a  CH  in  the  Unified  Soil  Classification 

System.  Itom  a  dqpth  of  27  to  32  ft,  an  interbedded  stratum  of  coarse 

sand  and  gravel,  which  was  encountered  in  boring  A,  is  shown  in  fig.  3i 

this  stratum  was  encountered  in  all  three  borings.  A,  B,  and  C,  but  at 

different  levels.  The  grain-size  distribution  of  a  sand  sample  from  thi.s 

layer  is  shown  in  fig.  2.  From  a  depth  of  32  ft  down  to  at  least  a  depth 

of  67  ft  the  soil  is  a  graj^  silty  clay.  Hie  watei*  content  in  this  layer  is 

about  23  percent,  and  the  unoonfined  compressive  strength  averages  about 

p 

1  ton/ft  .  Hie  liquid  limit  is  37  and  the  plastic  limit  is  I6;  therefore, 
this  clay  is  a  CL  j.n  the  Unified  Soil  Classification  l^stmn.  Several 
consolidation  tests  were  run  on  samples  from  this  stratum  and  the  analyses 
of  the  e  log  P  curves  indicate  that  the  layer  is  essentially  normally 
loaded  with  a  compression  index  of  0.25 . 


Conditions  in  the  Siqierseismlc  Region 

8.  In  the  superseismic  region  the  velocity  of  the  air  shock,  U  ,  as 
shown  in  fig.  4,  is  much  faster  themi  the  seismic  velocity  c  in  the  soil, 
such  that  at  relatively  shallow  depths  the  angle  6  which  the  front  of  the 
stress  wave  in  the  soil  makes  with  the  horizontal  is  relatively  small.  For 
purposes  of  analysis,  the  stress  wave  is  assumed  to  proi>agate  vertically 
down  the  soil  columm  illustrated  in  fig.  4  such  that  vertical  compression 
occurs  in  the  soil  with  essentially  zero,  or  very  limited,  lateral  strain. 
The  significant  property  of  the  soil  under  these  conditions  which  controls 
the  deformation  under  a  given  loading  is  the  dynamic  constrained  modulus 
(secant)  of  deformation  .  In  a  constrained  or  one-dimensional  condi¬ 
tion,  the  general  shape  of  the  stress-strain  curve  is  as  shown  in  fig.  5* 

At  low  stresses  the  stress-strain  curve  is  concave  downward,  but  as  the 
pressure  level  increases  the  curve  passes  through  a  point  of  inflection 
and  becomes  concave  upward.  The  stress  level  at  which  the  point  of  inflec¬ 
tion  occurs  is  a  function  of  the  natural  cementation  of  the  soil,  the 
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(b)  (d) 

Fig.  4.  laical  conditions  in  the  superseismic  region 


Fig.  5.  laical  stress-strain  relation  for  an 
undisturbed  soil  in  one-dimensional  conipression 


maximum  preconsolidatlon  load,  and  the  strain  rate  at  vhich  the  test  was 
conducted.  It  has  been  shown^  that  for  nearly  identical  samples  of  silt, 
the  point  of  inflection  in  a  dynamic  test  occurs  at  a  stress  of  nearly 
two  times  the  stress  at  which  the  inflection  point  occurs  in  a  static  test. 
It  is  quite  apparent,  then,  from  fig.  5  that  the  constrained  modulus  is 
very  much  a  function  of  the  pressure  level.  At  low  stresses,  say  less  than 
1  psi,  the  constrained  modulus  is  related  to  the  initial  tangent 
shown  for  the  stress-strain  cttrve  in  fig.  5>  and  could  be  estimated  from 
the  seismic  velocity  c  in  the  field  means  of  equation  1  where  p  is 
the  mass  density  of  the  soil. 

=  pc^  (1) 

As  the  pressure  level  increases,  however,  the  modulus  decreases,  until  some 
point  beyond  the  point  of  inflection  where  the  constrained  secant  modulus 
begins  to  increase  with  pressure.  It  is  imperative,  therefore  (since  the 
modulus  is  so  dependent  upon  the  stress  level),  that  the  procedure  used 
for  determining  the  constrained  modulus  of  undisturbed  samples  from  n  given 
site  be  conducted  at  a  pressure  level  commensurate  with  the  pressure  level 
expected  in  the  field  at  that  location.  The  dynamic  testing  program  on 
the  WES  saiiq)les  was  planned  so  that  the  stress  levels  employed  would  sup- 

■i 

plement  the  test  results  reported  by  Davisson  and  Maynard"^  in  order  that 
constrained  modulus  data  would  then  be  available  on  the  Suffield  Experi¬ 
mental  Station  soils  over  a  wide  range  of  pressure  levels. 

Results  of  Dynamic,  Constrained  Modulus  (Secant)  Tests 

9.  The  Wilson  vibration  apparatus  was  utilized  to  test  16  samples 
of  undisturbed  soil  from  the  Suffield  Experimental  Station  from  depths 
ranging  from  2  to  60  ft.  Ihe  results  of  these  tests  were  reported  by 
Shannon  and  Wilson  and  are  shown  in  fig.  6.  This  apparatus  essentially 
gives  a  modulus  of  the  soil  material  at  very  low  stress  increments;  as 
shown  in  fig.  6,  the  modulus  as  determined  by  this  type  of  test  varies 
from  about  24,000  psi  at  the  surface  to  13,000  psi  at  a  depth  of  10  ft, 
and  from  a  depth  of  10  to  ^9  ft  is  nearly  constant  between  the  values 
of  12,000  and  1^,000  psi. 
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constrained  modull-depth  relations 


10.  The  constrained  dynardc  modtili  as  back-calculated  from  the 
seismic  velocities  by  equation  1  are  also  shown  in  fig.  6.  From  0  to  23  ft 
the  seismic  velocities  indicate  a  constrained  modulus  of  24,000  psi.  From 
a  depth  of  23  to  70  ft  the  seismic  velocities  indicate  a  constrained 
modulus  which  ranges  from  60C/,000  to  900 >000  psi.  The  discontinuity  in 
seismic  velocities  at  a  depth  of  23  ft  indicates  that  at  the  time  of  the 
seismic  survey  the  water  table  was  at  a  dq)th  of  23  ft.  It  should  be  noted 
that  the  constrained  moduli,  as  back- calculated  from  the  seismic  values  sure 
somewhat  higher  than  the  constrained  moduli  as  given  by  the  Wilson  vibra¬ 
tion  tests  because  the  stress  levels  involved  in  the  seismic  pulse  are 
even  lower  than  the  stress  levels  utilized  in  the  Wilson  vibration  tests. 

11.  Eynamic,  confined  compj.*e8sion  tests  were  also  conducted  at  the 
University  of  Illinois  for  the  Air  Force  Weapons  Laboratory."^  The  results 
of  these  tests  are  also  shown  in  fig.  6.  The  dynamic,  confined  compression 
tests  were  conducted  with  an  apparatus  as  shown  in  fig.  7;  in  these  tests 
the  sample  is  confined  by  a  steel  ring  and  loaded  dynamically  with  a  gas 
loader  which  subjects  the  sample  to  stress  levels  comnensurate  with  those 
expected  in  the  field.  Idle  dynamic  constrained  moduli  as  determined  from 
these  tests  are  shown  in  fig.  6  for  the  200- ,  1^0-,  and  100-psi  stress 
levels.  At  an  overstress  of  from  200  to  300  psi,  the  modulus  as  determined 
from  the  dynamic,  confin<»d  compression  tests  is  about  3000  psi  for  a  depth 
of  0  to  13  ft.  Idle  constrained  dynamic  modulus  at  130-psi  overpressure 
Increases  with  depth,  as  shown  in  fig.  6,  from  a  value  of  4300  psi  at  a 
depth  of  3  ft  to  a  value  of  12,000  psi  at  13  ft.  At  an  overpressure  of 
100  psi,  the  constrained  modulus  also  increases  with  depth  from  a  value  of 
8000  psi  at  a  depth  of  3  ft,  10,000  psi  at  a  depth  of  6  ft,  13 >000  psi  at 

a  depth  of  9  to  24,000  psi  at  a  depth  of  13  ft.  Idiese  data  Indicate 
clearly  that  as  the  stress  level  Increases  at  a  given  depth  the  constrained 
modulus  decreases;  and  that  for  any  given  value  of  the  vertical  stress,  o^  > 
the  constrained  modulus  generally  increases  with  depth. 

12.  Idle  constrained  modulus  as  calculated  from  the  propagation  veloo 
Ity  of  the  peak  stress  in  a  wave  propagation  test  conducted  at  the  VSS*  is 

*  Feasibility  of  aiimll  Blast  Load  Generator  Wave  1  rqpagatlon  Velocity 
Measurements  on  Undisturbed  Soil  Samples,  mBmorandum  for  Chief,  Nuclear 
Weapons  Effects  Division,  U.  AziQr  linear  Waterways  Experiment 
Station,  Vicksburg,  Mississippi,  March  1963> 
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also  shown  in  fig.  6.  This  test  was  conducted  on  a  sangjle  from  a  depth  of 
2  ft  at  a  stress  level  of  250  psi;  the  results  Indicated  an  effective 
modulus  of  2900  psi  which  is  in  good  agreement  with  the  constrained  modu¬ 
lus  values  at  this  stress  level  as  determined  from  the  dynamic^  confined 
conipression  tests  described  by  Davisson  and  Maynard.  The  general  objec¬ 
tive  of  this  test  was  to  measure  the  peak  stress  wave  velocity  and  residual 
strain  in  an  SBLG  test  on  a  5-in. -diameter  undisturbed  sample  utilizing  a 
dynamic  pressure  pulse  of  the  same  magnitude  as  encountered  in  the  field. 
The  undisturbed  sample  was  placed  in  a  sonicircular  cradle,  and  the  card¬ 
board  and  wax  protection  were  partially  trimmed  away,  as  shown  in  fig.  8a, 
to  form  a  "stacked  ring  device”  suround  the  saaple  in  order  to  minimize 
friction  losses  during  the  dynamic  loading  in  the  SBLG.  At  locations  shown 
in  fig.  8b,  Illinois  Institute  of  Technology  Research  Institute  (lITRl), 
miniature  stress  gages  were  inserted  into  cavitites  carved  in  the  center 
of  the  sanple  as  shown.  The  sample  was  then  placed  in  the  SBLG  and  sur¬ 
rounded  with  two  sheets  of  teflon,  which  were  in  turn  surrounded  by  dense 
sand.  During  placement  of  the  sanple  in  the  generator,  IITRI  pressure 
cells  were  placed  at  the  top  and  bottom  of  the  saaple  and  a  Road  Resecurch 
Laboratory  (RRL)  cell  was  placed  in  the  sand  at  a  level  corresponding  to 
the  base  of  the  saaple.  The  bonnet  of  the  SBLG  was  then  secured  over  the 
specimen,  and  the  generator  was  fired  at  a  pressure  of  250  psi.  Hie  test 
results  indicate  that  the  average  velocity  of  the  peak  stress  in  the  top 
0.86  ft  of  the  Snowball  specimen  was  370  fps.  (Hie  peak  stress  velocity' 
in  the  sand  was  1080  fps. }  The  peak  stress  wave  velocity  was  laeasured 
in  the  Snowball  exp^iment^  from  times  of  arrival  on  pressure  gages  ar¬ 
ranged  in  a  vertical  column  near  the  BOO-psl  range  (boring  C).  Hie  field 
measurements  Indicated  that  the  peak  stress  velocity  was  390  fps  in  the 
iQper  13  ft.  Hius .  the  VES  wave  propagation  test  on  an  undisturbed  sample 
yielded  a  peak  6..ress  wave  velocity  for  the  Suffield  soil  which  was  in 
good  agreement  with  field  observations.  Also  the  effective  modulus  back- 
calculated  from  the  velocity  of  the  peak  stress  is  in  good  agreement  with 
the  constrained  modulus  determined  at  comparable  pressure  levels  in  the 
dynamic,  confined  coapressico  tests  reported  by  Davisson  and  Maynard.^  It 
is  therefore  concluded  that  the  type  of  wave  poropagation  test  ccmducted 
at  WES  looks  quite  promising  and  should  be  investigated  further. 
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13-  The  aforementioned  test  results  and  discussions  demonstrate  the 
dependence  of  the  effective  constrained  modulus  on  the  pr-'ssure  level. 

Both  the  field  observations  and  laboratory  tests  conducted  at  higher  over¬ 
pressures  demonstrate  the  stress  vave  velocity  for  a  peak  stress  between 
200  apd  300  psi  to  be  about  k/lO  of  the  seismic  wave  velocity  for  the  top 
13  ft  of  this  soil  profile.  These  data  then  only  serve  to  demonstrate  that 
the  seismic  velocity  is  at  best  only  an  index  property  of  the  profile  which 
serves  to  give  em  tipper  limit  to  the  velocity  of  propagation  and  em  t^tper 
limit  to  an  effective  constrained  modulus. 

lU.  Ihese  data  indicate  that  constrained  moduli  values  back- 
calculated  from  the  seismic  velocity  could  be  in  error  by  as  much  as  a 
factor  of  8  for  those  shallow  portions  of  the  soil  profile  subjected  to 
an  overpressure  of  200  to  300  psi.  An  exaaple  follows  which  dePKmstrates 
the  use  of  a  procedure  for  selecting  the  constrained  modulus  with  depth 
at  a  given  location  if  the  peak  side-on  overpressure  at  the  surface  can 
be  estimated  prior  to  the  shot. 

Selection  of  a  Constrained  Modulus-Depth  Relation 


15.  Since  the  effective  constrained  modulus  at  a  given  depth  is  a 
function  of  the  peak  vertical  stress  at  that  depth,  the  first  step  in  the 
selection  of  a  constrained  modulus-depth  relation  is  to  estimate  the 
attenuation  of  vertical  stress  with  depth.  Predicting  the  attenuated  peak 
vertical  stress  with  depth  has  been  the  object  of  considerable  study  and 
will  not  be  discussed  in  detail  here;  the  method  used  will  be  essentially 
that  given  in  the  Air  Force  Design  Nsnual.*^  It  Is  well  established  that 
this  method  gives  results  which  sre  In  fairly  close  agreement  to  pressures 

D 

observed  In  the  field  even  though  a  recent  study  has  shown  the  assuaptlons 
underlying  this  method  to  be  questionable. 

16.  In  this  method,  the  peak  vertical  stress  at  depth  a  can 
be  expressed  as  a  function  of  the  peek  slde-oo  overpressure  and  the 
attenumtloo  factor  as  given  In  equation  2 
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(3) 


Where,  a  is  given  by  equation  3 
z 

-  1 

“z  "  1  +  z/l 
'  w 

and  is  given  by  equation  U 

=  .30  n  (^f 

BO  '  ' 

17-  A  curve  of  the  calculated  peak  attenuated  vertical  stress  with 
depth  is  shown  in  fig.  9b  for  a  location  200  ft  from  GZ,  which  corresponds 
to  WES  borlxig  C.  At  this  range  the  predicted  value  of  the  peak  slde-on 

9 

overpressure  by  was  310  psl.  In  fig.  9^  are  shown  the  various  data 
which  define  the  variation  of  the  dynamic  constrained  modulus  with  depth 
for  the  Suffleld  Experimental  Station  soils  as  determined  by  the  test 
procedures  described  previously  which  utilize  various  levels  of  stress. 

Ihe  task  then  is  to  look  at  the  data  available  on  the  soil  in  fig.  9^  and 
the  predicted  attenuatea  vertical  stress  in  fig.  9h  and  construct  a  con¬ 
strained  modulus -depth  curve  ^diich  accounts  for  the  pressure  level  expected 
at  a  given  depth.  For  exaiqple,  the  dynamic  test  results  have  shown  that 
saoples  taken  near  the  surface  and  subjected  to  dynamic  overpressures  in 
test  devices  from  200  to  300  psi  will  have  a  constrained  modulus  of  about 
3000  psl.  Therefore,  at  the  200-ft  station,  which  has  a  peak  overpressure 
at  the  surface  of  310  psi,  elements  of  soil  near  the  surface  will  be  ex¬ 
pected  to  have  a  constrained  modulus  of  3000  psi.  At  a  depth  of  10  ft 
at  the  200-ft  station,  the  peak  attenuated  stress  as  given  by  fi8>  9h  la 
approximately  130  psi.  Therefore,  at  a  depth  of  10  ft  on  fig.  9*  the 
constrained  modulus  of  interest  would  be  7000  psi ,  the  value  of  the  con¬ 
strained  modulus  where  the  130-psl  line  passes  through  a  depth  of  10  ft. 

At  deeper  depths  (13  to  23  ft), the  constrained  modulus  should  approach 
the  1A,000-  to  13,000-psi  value  as  given  tgr  the  Shannon  and  Wilson  vibra¬ 
tion  test,  becau  it  has  been  shown  previously that  the  Wilson  vibration 
test  gives  a  reliable  value  of  the  constrained  modulus  when  the  strain  re¬ 
covery  approaches  100  percent.  Since  the  degree  of  saturation  in  the  13- 
to  23-ft  depth  rai^e  is  over  96  percent,  the  residual  strains  for  vertical 
stresses  approaching  100  psi,  as  is  the  case  at  this  station,  are  very 
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small.  Hie  constrained  modulus-depth  relation  arrived  at  hy  the  vrit*ir 
for  this  particular  situation  is  shown  by  the  smooth  black  line  in  fig.  9^^. 

This  relation  is  smooth  down  to  a  depth  of  23  ft  where  the  water  table 
table  would  cause  an  abnqpt  discontlnui'ty  and  the  constrained  modulus 
would  Junp  to  a  level  of  300^000  psi  or  more.  Since  the  stress  levels  ia  ' 

this  experiment  were  very  low  below  a  depth  of  23  ft,  when  considered  with 
a  constrained  modulus  of  300,000  or  more,  the  strains  below  a  depth  of 
23  ft  would  be  Insigniflcemt  in  causing  displacements  and  were  not  con¬ 
sidered  in  this  Einalysls. 

18.  Hie  peak,  attenuated,  vertical  stress-depth  curve  at  a  range  of 
250  ft,  which  corresponds  to  WES  boring  D,  is  shown  in  fig.  10b.  The  con¬ 
strained  modulus-depth  relation  is  given  by  the  smooth  black  curve  in 
fig.  10a  and  was  arrived  at  by  reasoning  and  Judgment  similar  to  that 
described  above. 

! 

Airblast-Induced  Ground  Motion  Predictions 

19.  The  peak  transient  vertical  displacements  at  the  surface  euad  a 
depth  of  5  ft  have  been  conputed  for  the  200-  and  250-ft  stations.  Hie 
calculations  are  based  vpon  the  predicted  peak  overpressure  at  these  ranges 
as  given  by  reference  9  and  the  oveipressure-time  curves  as  shown  in' 

fig.  11,  which  are  idealizations  of  data  given  by  BRL.*  The  peak  vertical 

stress-depth  relation  was  computed  from  equations  2,  3,  and  4  by  the  pro- 

7 

cedure  given  in  the  Air  Force  Design  Manual;  the  attenuated  vertical 
stress-depth  relations  are  shown  in  figs.  9  and  10. 

20.  As  the  airblast  arrives  at  a  point  on  the  surface  In  the  super- 
seismic  region,  an  alrblast-lnduced  stress  wave  begins  to  propagate  down¬ 
ward;  as  it  propagates  downward  the  peak  vertical  stress  attenuates  and 
the  rise  time  of  the  wave  Increases  with  depth  as  shown  in  fig.  ^  (page  6). 

The  rise  time  Increases  because  the  peak  stress  travels  at  a  velocity 
lower  than  the  seismic  velocity;  the  seismic  velocity  is  related  only  to 


*  Overpressure  Data  Obtained  by  Project  1.1  During  Operation  Snowball. 
conmunlcation  from  U.  S.  Anty  Ballistic  Research  Laborntories  to 
Director,  U.  S.  Amy  Engineer  Vateiways  Experiment  Station,  Vicksburg, 
Mississippi,  17  November  1964. 
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the  initial  surrival  of  the  pulse.  For  purposes  of  predicting  vertical 
ground  displacements  in  this  study  >  the  position  of  the  stress  wave  at 
times  ,  Tg  ,  and  ,  as  shown  qualitatively  in  fig.  Ud,  were  deter¬ 
mined.  At  any  given  time,  say  ,  the  variation  of  stress  with  depth 
is  assumed  to  be  linear  between  the  following  critical  points:  zero  stress 
at  the  depth  of  the  wave  front,  the  attenmted  peak  vertical  stress  at 
the  depth  of  the  peak  stress,  and  the  overpressure  at  the  siirface  for  time 
as  shown  in  fig .  4c .  At  various  times  the  strains  within  the  depth 
of  wave  propagation  are  integrated  to  give  the  displacement  at  the  surface. 
The  details  of  this  procedure  are  outlined  in  reference  11. 

21.  For  the  purpose  of  predicting  ground  motion  in  this  study,  only 
the  upper  23  ft  was  considered  because  the  soils  above  this  d^th  were 
rather  compressible  and  below  a  depth  of  23  ft,  which  corresponds  to  the 
depth  of  the  water  table,  the  soil  had  a  constrained  modulus  greater  than 
300,000  psij  therefore,  very  small  strains  would  result  below  this  dq>th 
for  the  pressure  levels  commenstirate  with  the  200-  to  250-ft  range.  Qhe 
initial  pulse  was  assumed  to  propagate  at  a  velocity  of  1000  fps  and  the 
peak  stress  was  assumed  to  propagate  at  ^00  fps  iu  the  upper  23  ft.  After 
the  vertical  stress-depth  relations  were  determined  at  various  times,  as 
shown  q\ialltatlvely  in  fig.  4d,  the  constrained  modulus-depth  relations 
from  figs.  9  and  10  were  used  with  the  vertical  stress-depth  relations  to 
arrive  at  a  strain  distribution  with  depth  at  various  times.  Ihe  strains 
were  then  Integrated  to  give  the  displacements  of  convenient  layers,  as 
indicated  by  layers  A,  B,  C,  D,  etc.,  on  the  right  side  of  figs.  9  and  10. 
By  summing  the  displacements  of  the  layers  below  the  surface  and  below  a 
depth  of  5  ft,  the  displacements  at  the  sivface  and  at  a  depth  of  ^  ft 
were  determined.  !Ihe  results  of  these  preliminary  calculations  yield  upper 
bound,  peak,  transient  surface  displacements  at  the  200-  and  250-ft  sta¬ 
tions  of  8.0  +  0.5  in.  and  5*4  +  0.5  in.,  respectively,  >rtille  the  calcu¬ 
lated  peak,  transient  displacements  at  a  depth  of  5  ft  at  the  same  two 

stations  are  4.5  and  2.2  in.,  respectively.  Ground-motion  measurements 
12 

by  Murrell  indicate  that  the  peak,  transient  displacements  at  a  depth  of 
5  ^/were  4.9  and  2.1  in.  for  the  200-  and  250-ft  stations,  respectively, 
which  are  in  reasonable  agreement  with  the  c<..  ted  values. 
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Conclttalons 

22.  The  dynamic  constrained  moduli  data  available  on  the  Suf field 
Experimental  Station  soils  definitely  indicate  that  it  is  of  prime  impor¬ 
tance  to  measure  uhe  constrained  modulus  of  a  soil  at  the  stress  levels 
expected  at  the  location  of  interest  in  the  field.  Direct  inference  of 
the  constrained  modulus  from  the  seismic  velocity  or  vibration  tests  will 
greatly  overestimate  the  constrained  modulus  and  result  in  an  underestima¬ 
tion  of  the  ground  motion. 

23.  Ihe  wave  propagation  test  conducted  at  WES  on  an  undisturbed 
sample  of  Suffleld  soil  gave  the  same  velocity  of  the  peak  stress  as 
observed  at  the  same  depth  in  the  field;  the  constrained  modulus  back- 
calculated  from  this  peak  stress  velocity  agrees  quite  well  with  the 
dynamic  constrained  modulus  determined  at  the  same  stress  level  in  dynamic 
confined  compression  tests  on  the  same  soil  at  the  University  of  Illinois. 

24.  The  ground  displacements  predicted  on  the  basis  of  the  dynamic 
properties  of  the  soil  presented  herein  agree  reasonably  well  with  the 
observed  values.  It  is  emphasized,  however,  that  the  analysis  presented 
is  preliminary  and  further  analysis  of  the  observed  ground  notions  in 
terns  of  the  dynamic  soil  properties  is  needed.  Additional  analyses 

of  the  ground-motion  data  will  be  made  in  the  final  report  on  the 
earth  motion  study  (in  publication). 
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alal  rapart  atedbar  by  aiUab  *a  daaoMai  ■W  ba  Mn 
aid  caaitaMad  by  iba  arttlaailab  actlvtoy.  TblanaMI 
ba  aatbaa  to  tWa  lapan. 

«b  OTimiWPOlirimaejICEb  VlbarBanbM 
iMMlNMd  any  «tbar  tepan  aaabara  faMter  bf  tea  arm 
or  by  dto  cpaaaarA  alsa  atear  Uda  a«bbai(ab 


to.  AVAlLABUTY/LDirrATIONIIOTtCBS;  Ealar  aay  Iba- 
Itailoaa  aa  fotbar  dlaaaataailoa  af  tba  lapart,  otbar  thaa  tboaa 
iayoaad  by  aacarity  clacaiOcatton,  aatet  alaadaid  atatoateota 
ouch  at: 

(1)  “Owllflad  roaaaalua  cMy  ablata  coplao  af  tbia 
rapoit  bum  DDC'* 

(2)  "Faralgt  iaanaacMite  aad  diaaaalaatiaa  of  ttea 
lopait  by  DDC  la  aai  aalbailaaA** 

(3)  “U  S.  ftnvataatete  cpaaclaa  aay  ablala  caplan  of 
Ibla  lapan  dbocUy  baa  MIC.  Otbar  paallSad  DDC 
uaara  abcll  lapraat  tbraapb 


<4)  “U.  E  aUnary  apaadoa  aay  obtain  coplaa  af  tela 
rapart  dbactly  baa  DOC  Otea  ptellWad  aaara 


<$)  **  AU  dteMbaitaa  af  Uda  rapart  la  aaadtellaE  Qaal* 
iflad  DDC  naaa  abaU  laaaaai  ttaaaab 


If  tea  rapart  baa  baaa  ftoafabtd  to  tba  OWteaafTaabateal 
Sarvieaa,  Dip  atonal  af  CanaacA  lar  aaia  la  tea  ptelte.  ladb 
caa  Uda  fact  aad  arda  tea  prica  If  baaaa 

IL  SUPPLOIBirTARV  NOTBk  Uaa  far  addltlaaU  aiplaaa 


iPOlWOmiWO  imUTiy  ACnyiTTi  laartB 
:  dapartaaalal  pcalaat  afEaa  ar  libaralaiy  teaaaarl 
I  far)  tba  tiaaacb  aad  davalapaatb  hwtada  adtea 


IX  ABEntACT  Bidnr  aa  cbawact  tl*l*S  a  brief  aad  factaal 
aaaaiiy  af  tea  dacaaaal  ladicattva  af  Iba  laaart.  avaa  teaatb 
It  aay  alaa  appaa  alaaadara  la  tea  bady  af  Hm  taabrUeal  ra* 

b#  —  ^  k  —  —  —  b  ia  ^ 

II  SpBCV  Iv  w 

aball  ba  altacbad. 

h  to  bbUdy  daalrablo  teal  tea  abatract  of  alaaalflad  ia> 
porta  ba  aacliittflad.  Saab  paraprapb  af  tea  ibabact  aball 
and  trite  an  (adicatlaa  af  tea  allltary  aacartly  elaaalltoMlao 
at  Iba  lafarrmtlta  to  tba  paravapb,  iipiaaaidid  aa  fTEX  fEX 

fCX  arflf). 

Tbara  la  an  Itodiailao  aa  tea  laateb  at  tea  abatract.  Harr- 
avat,  tec  anppaarad  laapib  la  baa  IN  la  211  aarda. 

lA  EET  WORDEt  Ray  aarda  arc  tarbatrally  aaaalatfal  taraa 


